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Point contacts provide simple connections between macroscopic particle reservoirs. In electric cir¬ 
cuits, strong links between metals, semiconductors or superconductors have applications for funda¬ 
mental condensed-matter physics as well as quantum information processing. However for complex, 
strongly correlated materials, links have been largely restricted to weak tunnel junctions. Here we 
study resonantly interacting Fermi gases connected by a tunable, ballistic quantum point contact, 
finding a non-linear current-bias relation. At low temperature, our observations agree quantitatively 
with a theoretical model in which the current originates from multiple Andreev reflections. In a wide 
contact geometry, the competition between superfluidity and thermally activated transport leads 
to a conductance minimum. Our system offers a controllable platform for the study of mesoscopic 
devices based on strongly interacting matter. 

PACS numbers: 


The effect of strong interactions between the con¬ 
stituents of a quantum nrany-body system is at the ori¬ 
gin of several challenging questions in physics. Whilst 
the ground states of strongly interacting systems are in¬ 
creasingly better understood JT], the properties out of 
equilibrium and at finite temperature often remain puz¬ 
zling, as these are determined by the excitations above 
the ground state. In laboratory experiments, strongly in¬ 
teracting systems are found in certain materials, as well 
as in quantum fluids and gases [lj. In solid-state sys¬ 
tems, a conceptually simple and clean approach to probe 
non-equilibrium physics is provided by transport mea¬ 
surements through the well-defined geometry of a quan¬ 
tum point contact (QPC) j2H3]. Yet, the technical hur¬ 
dles to realise a controlled QPC between strongly corre¬ 
lated materials pose a big challenge. Ultra-cold atomic 
Fermi gases in the vicinity of a Feshbach resonance, the 
so-called unitary regime, provide an alternative route to 
study correlated systems [5]- Superfluidity has been es¬ 
tablished at low temperature [B], but the finite- temper¬ 
ature properties are only partially understood HH9], a 
situation similar to the field of strongly correlated mate¬ 
rials. 

Recent progresses in the manipulation of cold atomic 
gases have allowed to create a mesoscopic device featur¬ 
ing quantised conductance between two reservoirs in the 
non-interacting regime [16] , We use this technique to 
create a QPC in a strongly interacting Fermi gas consist¬ 
ing of 1-7(2) x 10 5 6 Li atoms in each of the two lowest 
hyperfine states, in a magnetic field of 832 G, where the 
interaction strength diverges due to a broad Feshbach res¬ 
onance. The atoms form a strongly correlated superfluid, 
with a pairing gap larger than the chemical potential [3] . 
Typical temperatures in the cloud are T = 100(4) nK at 
a chemical potential of q = 360 nK • kp- The setup is 
presented in Figure Eh m- The QPC is characterised 


by transverse trapping frequencies of v x = 10.0(4) and 
v z = 10(3) kHz in x- and ^-direction. An optical at¬ 
tractive ’’gate” potential is used to tune the chemical 
potential and the number of channels in the QPC (see 
Figure 03) m- We prepare two atomic clouds (reser¬ 
voirs) with an atom number difference AN while blocking 
transport through the QPC with a repulsive laser beam. 
This results in a chemical potential bias between the two 
reservoirs A fi = f(AN,T/Tp), with the Fermi tempera¬ 
ture Tp and the function / derived from the equation of 
state (EoS) of the trapped Fermi gas at unitarity [8) ITT]. 

We open the QPC and measure AN as a function of 
time t. Figure presents this evolution for various 
strengths of the gate potential Vq- We observe that the 
AN decays from its initial value to zero over a time scale 
of 0.5 to 1.5s. The shape of the decay curves deviates 
from an exponential and is a direct manifestation of a 
non- linear relation between AN and the atom current. 

We extract the numerical derivative of these data HU. 
yielding the instantaneous current at a certain Ap, shown 
in Figure [2j3. We normalise the current and bias by the 
strength of the pairing gap A, using the known relation of 
A with the chemical potential /i for the low temperature 
Fermi gas, A = q/S, ■ q, with q = 0.44 [T^j and £ = 0.37 

mm- 

The current-bias characteristics in Figure [2j3 are 
strongly non- linear for all the choices of Vq, featuring a 
very strong response at low bias. The lrigh-bias regime 
approaches a linear dependance with a nonzero intercept 
on the current axis, marking an excess current that de¬ 
pends on the strength of the Vq- The current is much 
larger than what is observed for non-interacting atoms 
m, as observed in earlier measurements on strongly in¬ 
teracting atoms in multimode channels m- 

We model the experimental system as two superfluid 
reservoirs connected by a single particle hopping mech- 
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FIG. 1: Concept of the experiment. (A) Schematics of the two atom reservoirs (blue) connected by a QPC. Two repulsive 
beam (green) confine the center of the cloud in x- and 2 -direction, the attractive gate beam (red) tunes the density in the 
QPC. The dark contours are schematics of the beam profiles. (B) Potential landscape in the plane 2 = 0. Close to the QPC 
the attractive gate creates areas of high density (blue). (C) Theoretical model for the QPC. Both sites of the QPC have 
a defined atom number, imposing a chemical potential /xl, /xr and a pairing gap Al,Ar. The transparency a is an energy 
dependent function describing the transmission of single particles from one site to the other. (D) Transport via multiple 
Andreev reflections. Coherent tunnelling of pairs allows for the creation and tunnelling of a single particle excitation (pair 
breaking) leading to a DC current, even for A/x <C A. 



anism solely characterised by the transparency a of the 
QPC (see Figure [2p) 0HMB]. This model excludes 
any finite size and geometry-dependent effects, which we 
think of as absorbed in a. It is motivated by the large 
proximity effect in superfluids separated by a ballistic 
normal barrier |20l . Indeed, we expect a coherence length 
of ~ 3/mi, where Vp is the Fermi velocity, fcs is 
Boltzmann’s constant and T is the temperature of the 
gas. This is comparable to half the length of the chan¬ 
nel (5.6 /xm), thus we approximate the channel with a 
point-like connection. Using a non-equilibrium Keldysh 
Green function technique m with mean-field approxi¬ 
mation m we compare theoretical current-bias curves 
with our data. 

Since the pairing gap and the EoS of the unitary gas 
are known a priori, the only free parameter in our model 
is the transparency a n for each transverse mode n in 
the QPC HU. The solid lines in Figure [2j3 show the re¬ 
sults with the best fits of a. For the two lowest gate po¬ 
tentials we obtain good agreement with a single channel 
model, whereas for higher gate potentials three channels 
are required, in agreement with our reference measure¬ 
ment with a weakly interacting Fermi gas HU. 

The agreement between theory and experiment clar¬ 
ifies the microscopic origin of the current. Reflecting 
the strongly interacting nature of the system, we have 
A/x <C A. In this regime, a current flow is allowed via 
multiple coherent reflections of quasiparticles between su¬ 
perconducting reservoirs, i.e., multiple Andreev reflec¬ 
tions, illustrated in Figure [lp HS- The gap for a single 
particle transfer can be bridged by the simultaneous, co¬ 
herent transfer of n pairs if 2nA/x > A, with a probability 
of order a 2n . As is seen in Figure |2p, the drop of cur¬ 
rent observed at low bias corresponds to A/x ~ A(1 — a) 
m, where the finite transparency a suppresses the corre¬ 
sponding Andreev processes. In the very low bias regime, 
not resolved in the experiment, the DC current is expo¬ 


nentially suppressed and an oscillating current caused by 
the energy mismatch between the two reservoirs occurs. 
This averages to zero in the DC limit, and represents an 
AC Josephson current adding to the DC response HHI221- 

m- 

We now investigate the current-bias relation as a func¬ 
tion of temperature, for a fixed gate potential Vq = 
674 nK. To this end, we introduce a controlled heat¬ 
ing of both reservoirs before the transport is started, 
using variable amplitude parametric heating. With this 
method we explore a temperature range of 124-290 nK, 
from a deeply superfluid regime up to the superfluid-to- 
normal transition point. We measure the decay of parti¬ 
cle imbalance with increasing temperature and observe a 
crossover towards exponential decay when temperature is 
above 145 nK. We extract the current- bias characteristic 
(Figure |3]4) using the known finite-temperature EoS of 
the unitary Fermi gas emu. With increasing temper¬ 
ature the non-linearity disappears and the current glob¬ 
ally decreases. We interpret this as the disappearance of 
the superfluid contribution to transport as temperature 
is raised. 

From these data, the differential resistance R at low 
bias is estimated by fitting a line to the low bias region 
of the curve. The result is presented in Figure [3)3, where 
the decrease in R is clearly visible as temperature is de¬ 
creased. We compare this to a model independent mea¬ 
sure of the non-linearity of the characteristic provided 
by the y 2 parameter of an exponential fit to the entire 
decay curves HU- This parameter, as well as the fitted 
timescale r of the exponential decay tracks the measured 
differential resistance (Figure [3p). 

While in the low temperature regime, non-linearity is 
captured by our mean- field model, the model fails to re¬ 
produce the resistance at high temperature, indicating a 
breakdown of the mean-field description. In particular, 
it predicts a resistance in the linear regime given by the 
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FIG. 2: Non-linear characteristic of the QPC. (A) Time evolution of the particle imbalance A N/N for V g =341 nK (filled 
circle), 443 nK (open square), 544 nK (filled triangle), 645 nK (open diamond) and 747 nK (filled hexagon). (B) Current-bias 
characteristics normalised with respect to A. The error bars represent the variation of three averaged data sets. Negative 
values of the current are artifacts from the numerical derivation process. The red lines show the result of Keldysh calculations 
with the transparency a of the QPC as the only free parameter (see ED- For clarity the curves are shifted vertically by 4 
units. 


non-interacting Landauer formula, while the resistance 
that we measure is one order of magnitude smaller. This 
constitutes an indirect evidence that the high tempera¬ 
ture state of the gas is not a Fermi liquid, since Fermi 
liquid theory leads to the Landauer formula, setting the 
upper limit for the current carried by an ideal contact. 
One possible explanation is the presence of superfluid 
fluctuations, which are expected to be large close to Tq 
[2S1- Indeed, it is known that a one-dimensional Fermi 
gas held between attractively interacting leads, a proto¬ 
type of non-Fermi liquid system [29j , shows an enhanced 
conductance 

The physical picture emerging from these measure¬ 


ments relies on the finite- temperature properties of the 
reservoirs. Some of the intrinsic physics of the channel 
appears when the mode spacing in the QPC is compa¬ 
rable to the temperature range explored. In this situa¬ 
tion, several modes are thermally populated [33], enhanc¬ 
ing transport for increasing temperature and competing 
with the subsequent reduction of the superfluid current. 
To show this, we decrease the confinement in the QPC 
along the direction to 5 kHz and systematically mea¬ 
sure the current at the largest bias as a function of tem¬ 
perature for various gate potential. The results are shown 
in Figure [4]4 where the current at high bias normalised 
to the bias is shown as a function of temperature for var- 
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FIG. 3: Finite-temperature transport properties. (A) Current -bias curves for temperatures 124 nK (open circle), 135 
nK (filled triangle), 145 nK (open triangle), 170 nK (filled hexagon), 223 nK (empty diamond) and 290 nK (filled square). 
For better visibility the curves have been shifted vertically by 40 units. Solid lines show best fits of the Keldysh calculations. 
Dashed lines are linear fits to the low bias regime, which yields the differential resistance R = shown in (B) as a function 

of decreasing temperature (C) Decay time r of the particle imbalance relaxation (open circle) obtained by fitting decay 

curves with an exponential; y 2 of the corresponding exponential fits (filled diamond) as a function of decreasing temperature 
The gray solid line marks the superfluid transition at k ^ T =3.56 




FIG. 4: Competition between superfluidity and thermally activated transport. (A) Normalised current as a function 
of temperature for various Vg■ 0 nK (open circle), 12 nK (filled triangle left), 54 nK (open hexagon), 97 nK (filled triangle up), 
139 nK (open diamond), 223 nK (filled triangle down), 308 nK (open pentagon), 519 nK (filled triangle right), 942 nK (open 
square). Inset: Zoom into the transient region for three selected values of Vg- (B) The same data as a function of The 
gray solid line marks the superfluid transition at k ^ T = 3.56. Few representative error bars are shown. 


ious gate potentials. In the linear response regime, this 
current-bias-ratio reduces to the differential conductance. 

For large Vg , we observe a decrease of the current 
with temperature over one order of magnitude. In this 
superfluid-donrinated regime, the pairing gap is large 
compared to both temperature and level spacing in the 
QPC, analogous to the previous measurements. At low 


Vg, the pairing gap is small and vanishes upon heating. 
In this channel-dominated regime, we observe an increase 
of the normalised current with temperature, which we at¬ 
tribute to thermal activation of transport channels in the 
QPC. There again, the current in the high temperature 
regime is much higher than predicted by a Fermi- liquid 
based Landauer formula. 
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At equilibrium in the reservoirs, superfluidity is uni¬ 
versally related to the local fugacity [34] due to the scale 
invariance of the unitary Fermi gas. This suggests the 
parameter pi/ksT as a common dimensionless scale for 
comparing conductances at various temperatures and Vq- 
Figure [4)3 presents the normalised current as a function 
of n/ksT. The data sets showing decreasing current with 
increasing temperature are all grouped in the high fu¬ 
gacity regime, below the expected superfluid transition 
point, confirming that this regime is dominated by su¬ 
perfluidity. Conversely, in the low fugacity regime the 
current increases with temperature, corresponding to the 
channel dominated regime. The crossover takes place 
close to the same fugacity for all the gate potentials, and 
is close to the universal transition point for the unitary 
Fermi gas at the center of the cloud. We expect that the 
exact location of the crossover, as well as the conduc¬ 
tance at the minimum depend on the details of the chan¬ 
nel geometry such as its energy dependent mode spacing. 
In addition, proximity effects should be reduced at high 
temperature and one dimensional physics could emerge in 
the QPC making the results dependent on its length [ 35 ]. 
Our setup, allowing for a direct and independent control 
of the geometry, could be used to investigate such effects 
in future experiments. 
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[1] A. J. Leggett, Quantum liquids: Bose condensation 
and Cooper pairing in condensed-matter systems (Oxford 
University Press, 2006). 

[2] N. van der Post, E. T. Peters, I. K. Yanson, J. M. van 
Ruitenbeek, Phys. Rev. Lett. 73, 2611 (1994). 

[3] E. Scheer, P. Joyez, D. Esteve, C. Urbina, M. H. Devoret, 
Phys. Rev. Lett. 78, 3535 (1997). 

[4] O. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod, 
C. Renner, Rev. Mod. Phys. 79, 353 (2007). 

[5] W. Zwerger, The BCS-BEC crossover and the unitary 
Fermi gas, vol. 836 (Springer Science & Business Media, 
2011 ). 


[6] M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C. H. 
Schunck, W. Ketterle, Nature 435, 1047 (2005). 

[7] S. Nascimbene, et al., Phys. Rev. Lett. 106 , 215303 

( 2011 ). 

[8] M. J. H. Ku, A. T. Sommer, L. W. Cheuk, M. W. Zwier¬ 
lein, Science 335 , 563 (2012). 

[9] Y. Sagi, T. E. Drake, R. Paudel, R. Chapurin, D. S. Jin, 
Phys. Rev. Lett. 114 , 075301 (2015). 

[10] S. Krinner, D. Stadler, D. Husmann, J.-P. Brantut, 
T. Esslinger, Nature 517 , 64 (2015). 

[11] Materials and methods are available as supplementary 
material on Science Online. 

[12] A. Schirotzek, Y.-i. Shin, C. H. Schunck, W. Ketterle, 
Phys. Rev. Lett. 101 , 140403 (2008). 

[13] G. Ziirn, et al, Phys. Rev. Lett. 110 , 135301 (2013). 

[14] D. Stadler, S. Krinner, J. Meineke, J.-P. Brantut, 
T. Esslinger, Nature 491 , 736 (2012). 

[15] G. E. Blonder, M. Tinkham, T. M. Klapwijk, Phys. Rev. 
B 25, 4515 (1982). 

[16] D. Averin, A. Bardas, Phys. Rev. Lett. 75 , 1831 (1995). 
[171 J. C. Cuevas, A. Martm-Rodero, A. L. Yeyati, Phys. Rev. 

B 54, 7366 (1996). 

[18] C. J. Bolech, T. Giamarchi, Phys. Rev. Lett. 92, 127001 
(2004). 

[19] C. J. Bolech, T. Giamarchi, Phys. Rev. B 71 , 024517 
(2005). 

[20] A. M. Zagoskin, Quantum theory of many-body systems 
(Springer, 1998). 

[21] A. Kamenev, Field theory of non-equilibrium systems 
(Cambridge University Press, 2011). 

[22] O. Avenel, E. Varoquaux, Phys. Rev. Lett. 60, 416 
(1988). 

[23] J. C. Davis, R. E. Packard, Rev. Mod. Phys. 74 , 741 

( 2002 ). 

[24] M. Albiez, et al, Phys. Rev. Lett. 95, 010402 (2005). 

[25] A. Ramanathan, et al., Phys. Rev. Lett. 106, 130401 

( 2011 ). 

[26] L. J. LeBlanc, et al, Phys. Rev. Lett. 106, 025302 (2011). 

[27] F. Jendrzejewski, et al, Phys. Rev. Lett. 113 , 045305 
(2014). 

[28] B. Liu, H. Zhai, S. Zhang, Phys. Rev. A 90, 051602 
(2014). 

[29] T. Giamarchi, Quantum physics in one dimension (Ox¬ 
ford University Press, 2003). 

[30] D. L. Maslov, M. Stone, Phys. Rev. B 52, R5539 (1995). 

[31] I. Safi, H. J. Schulz, Phys. Rev. B 52, R17040 (1995). 

[32] V. V. Ponomarenko, Phys. Rev. B 52, R8666 (1995). 

[33] J.-P. Brantut, et al., Science 342, 713 (2013). 

[34] T.-L. Ho, Phys. Rev. Lett. 92, 090402 (2004). 

[35] H. P. Birchler, V. B. Geshkenbein, G. Blatter, Phys. Rev. 
Lett. 92 , 067007 (2004). 



6 


Materials and Methods 
Theoretical model 


In general, transport properties of fermionic superfluid junctions with a few conduction channels can be dealt with 
in two different approaches: scattering formalism based on Bogoliubov-de Gennes equations [S1 j. !S2j and Hamiltonian 
formalism with nonequilibrium Keldysh Green function techniques jS3HS5j . We adopt the latter approach, which 
allows us to compute currents with high accuracy and potentially to treat more general systems with interaction and 
disorder. 

As discussed in the main text, since the junction is short, we approximate the ID channels connecting two reservoirs 
by a quantum point contact. To be specific, we consider the following Hamiltonian: 


H = H L + H R + H t , (1) 

Ht = - ^2 £ ^i jCT (0)tW(0) + h.c. (2) 

n <T=t,! 


Here, Hj (j = L or R) is the Hamiltonian of the reservoir j and is idential to the so-called single channel model 
)S6j . which consists of the quadratic part in the fermionic field rp including the kinetic energy, trapping potential, 
and chemical potential terms and the quartic part in rp describing the local pair interaction. On the other hand, 
H t describes tunneling processes between the reservoirs occurring at a single point denoted x = 0. The physical 
properties of the contact are assumed to be hidden in the tunneling amplitudes t n of each transverse mode of the 
ID contact n participating the tunneling. By assuming the above Hamiltonian and a steady state solution, one can 
obtain the current as follows: 


I = 


l, d(N R -N L ) , 
2 dt 


= £ 


itv 




L,cr 


(0)VtR,<x(0)} + h.C., 


(3) 


where N R and N R are the total number of particles in the reservoirs and the average (• ■ •) is expressed in terms of 
the Keldysh Green functions as discussed below. An important observation is the fact that the current depends only 
on the tunneling terms at position x = 0 and therefore one can obtain a local action by integrating out the position 
dependence of the reservoir Hamiltonians. 

While the above tunneling Hamiltonian has a simple structure, it is nevertheless difficult to obtain the current 
based on it since the Hamiltonian in each reservoir is quartic in 'll). To go further, we adopt a quadratic (mean-field) 
approximation, which renders the local Hamiltonian quadratic in ip as 

^ ^^ j ’ a ~ ( 4 ) 
rr ' ' 


Here, ji 3 and A j are the chemical potential and fermionic superfluid gap of the reservoir j defined locally. This 
approximation allows one to obtain the local retarded and advanced Green functions in each reservoir as follows 

EMU: 


r.a t \ 

9j M = 


— (u — Hi ± ir ]) 2 


— (w — Hj ± ir]) 

A, 


— (ui — Hj i iri) 


(5) 


where we introduce the Nambu representation for the fermions in each reservoir, ij = 0 + is an infinitesimal positive 
parameter that regularizes the Green functions, and the upper (lower) sign corresponds to the retarded (advanced) 
Green function, g r ( g a ). The Green functions are measured in units of an energy scale associated with the normal 
density of states at the Fermi level !S3l - iS5] . One obtains the Keldysh component of the Green functions by means 
of the fluctuation-dissipation theorem [ 57 ]: = (jjj — g’j) tanh((w — Hj)/ 2T) with the temperature T. By using the 

Green functions obtained above, the corresponding local action can be expressed as 


S = 



S t = 



£*n'0L, (7 MV’fl,<7M 

n,a 


+ h.c. 


( 6 ) 


(7) 
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Note that is a four components spinor consisting of the Keldysh and Nambu components, and the matrix g. 1 is 
given by 


9i 1 




( 8 ) 


where (g 1 -) 1 = -(fl£) 1 g k (g‘j) 1 EZi- 

The action © and the two reservoirs and the tunnelling term could be regrouped in an eight by eight matrix. It is 
in principle possible to invert it and thus to obtain the full solution of this problem. Note however that frequencies 
with equal positive and negative shifts from the Fermi level (±(w — n :j )) are used in the Nambu representation while 
absolute frequencies ui are conserved by the tunneling term (St). This implies that, even if it is quadratic, in the 
presence of a finite bias the action is no longer diagonal in frequencies and one must consider an infinite dimensional 
matrix to be inverted. However, only a discrete set of frequencies are coupled by the tunnelling term. Physically, 
such an infinite discrete set of frequencies represents multiple Andreev reflections (See Figure 1 D in the main text 
and Refs. [S1HS5] ). which are responsible for the non-linear I — V behavior in fermionic superfluids. Since the weight 
of the terms with frequencies far away from the Fermi levels decreases, it is possible to truncate the infinite series of 
the multiple Andreev reflections |S41 IS5| and thus to invert numerically the matrix with the desired accuracy. 

Inverting the matrix allows to compute the current. Indeed the current can be expressed as [S3HS5] 


I = E X / dujReGK ’ 


(9) 


where G K is the Keldysh component appearing in the inverse of the full action |S3HS5| . 

If both reservoirs are normal (A^ = 0), the offdiagonal elements in the Nambu representation disappear and the local 
action is diagonal in frequencies. One can then compute analytically the current and obtains the Landauer-Buttiker 

formula [551159] : 

J = E^ / du(fL(u>)-fR(w)], (10) 


where f :j is the Fermi distribution function in the reservoir j. and a n is the so-called transparency associated with 
the hopping amplitude t n as 





(ii) 


The transparency takes values in the interval [0,1], and perfectly transparent junctions corresponds to a = 1. 

Using the above formalism, one can thus compute the current for the nonzero A case as well by fixing the hopping 
amplitudes (see Figure [Slj). Results are given in Figure 2 B of the main text. In the case of the experimental system 
considered in the present paper, the junction is near perfect transparency awl. 


Data fitting 

The data presented in Figure 2B are fitted using the results of the numerical calculation of the current-bias charac¬ 
teristic, leaving only t n , the tunnelling amplitudes, as fit parameters. For each choice of gate potential the parameters 
are allowed to vary, reflecting the fact that the transparency is affected by the adiabaticity of the coupling and by the 
change of geometry induced by the gate potential. Figure [5T| presents the fitted tunnelling amplitudes as a function 
of the gate potential. The lowest fitted value is ~ 0.8, leading to a transparency for the corresponding mode larger 
than 0.95, compatible with the observation of quantized conductance in the non-interacting regime. Most of the fitted 
values are > 0.95, corresponding to a transparency a > 0.99. 


Thermodynamics 

We make use of the known thermodynamic relations for harmonically trapped and homogeneous unitary Fermi 
gases to extract the chemical potential bias and temperature of the gases. In this section we describe the employed 
procedure. 



Equation of states 


The density equation of state (EoS) of the homogeneous unitary Fermi gas can be expressed as a universal function 
of q = y/k B T [S10] 

n\ 3 = f„(q ) = —Li 3/2 (-e 9 ) F(q) (12) 


with the de Broglie wavelength A = y ■ Here the polylogarithm function represents the noninteracting part, while 

F(q) is the ratio n(q)/no(q) between the unitary and noninteracting density as measured in |S10j for —0.9 < q < 3.5. 
Using a fourth order virial expansion we extend the EoS to q < —0.9 IS101 ISli] . Further, we use the phonon model 
as described in [S121IS13] to expand the EoS to q > 3.5. The full EoS is then given by: 


fn(q) 


'Ei ><,j< '/ < -0.9 

F n (q) (—Li 3/2 (—e 9 )) , -0.9 < q < 3.5 


( 4?r ) 


3/2 
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5 / 2 ' 


3.5 < q 


(13) 


In the local density approximation (LDA), each point in the cloud has a local chemical potential given by 

q(x,y,z) = (mo~ (u 2 x x 2 + ut^y 2 + u) 2 z z 2 )^J (14) 

with the chemical potential yo in the bottom of the trap. The total atom number N in LDA is then given by the 
density integral over the full cloud 

N = J dxdydz f n ( q{x,y,z )) = ~^= M 0 (q 0 ) (15) 

with the average trapping frequency uj and the dimensionless moments |S14j 

/ Qo 

dq{qa~ q) (l+1)/2 f n (q) (16) 

-OO 

Using the known expression for the Fermi temperature ksTp = fuzifoN) 1 / 3 we obtain the degeneracy factor 

^=^M 0 (q 0 )^ 3 =:%o) (17) 

where the function h(q 0 ) formally represents the dependance of T/Tp on q 0 . We now invert h(q 0 ) and obtain the 
universal function g: 


- Mo J T \ 
Q0 k B T 



(18) 


Since we have determined T independantly, the function g directly yields the global chemical potential in the harmonic 
trap yo = g(T/Tp) ■ k B T in the reservoirs as a function of temperature. 

In a similar way we use the EoS to convert an atom number imbalance AN into a chemical potential bias Ay 
accounting for finite temperature. To this end we take the logarithm of Eq. (181 and derive it with respect to the 
atom number N 


AiogW = 

1 dy 1 9 1 (tf) t 1 
ydN ~ ~3 g Jjpj T^N 


(19) 

( 20 ) 


where we have denoted the derivative of g as g 1 . Using Eq. (18) to eliminate y and making a first 
in AN/N in Eq. (20), we obtain an expression for the chemical potential bias Ay as a function 


order approximation 
AN, N and T/T F 


Ay _ 1 , (T_\ 2 A N 

Kp ~ ~3 9 \T^J ~N~ 


( 21 ) 
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Finite temperature compressibility 


The finite temperature compressibility at unitarity is given by the derivative of Eq. (15): 

Comparing this to the zero-temperature limit compressibility 

ON 0 1 (6iV) 3/2 


C 0 = 


dpi 2 y^huj 


we find a ratio between C and Co of 


which is plotted in Fig. |S2| 


C VC M_ 2 (g(T/T F )) 

Co 12 2 / 3 7rV6 Mi ( ff ( t/T f )) 2/3 


( 22 ) 


(23) 


(24) 


Temperature calibration 


The temperature of the unitary Fermi gas can be obtained using the relation between the second moment of the 
density ( y 2 ) and the total energy per particle E/N |S15 


E = 3 mujy (j/ 2 ) 


(25) 


with the second moment 



dyni D (y)y 2 


(26) 


We select the direction y, because the cloud expands very slowly in y during the time of flight of 0.8 s, due to the 
magnetic confinement in y. The density profile nmiy) after time of flight is then indistinguishable from in situ 
pictures. Further, we use the EoS (Eq. ( fl2] )) to express the total energy per particle with reference to the Fermi level 
E/NEp. To this end we integrate the energy density e = n • (qo — q) over the full cloud and obtain 


E 


NEt. 


2tt\ 6 M 2 (q 0 ) 


^ Mt'\qo) 


(27) 


We can now combine Eq. (27) and Eq. (25) to obtain q 0 . Plugging this q 0 into Eq. (PL7|) yields the temperature T. 


Between the end of the transport process and the moment the cloud is imaged in the experiment cycle, the cloud 
is heated due to spontaneous emission and parasitic heating, with a heating rate R dependant on the power of the 
dipole trap laser. To obtain the temperature during transport T, we thus subtract the total heating in the dipole trap 
from the temperature at imaging Tj m . 


T = Tim — Rt (28) 

where t is the time span between the mean transport time and the imaging. In the data of Figure 2 the dipole trap was 
decompressed during transport and recompressed for imaging, which leads to a decreased heating rate. Furthermore 
the decompressed cloud has a Ep that differs from the one during imaging due to the dependancy on w. For an 
adiabatic recompression T/Tp stays constant, so the temperature scales with the ratio between the trap frequency 
of the compressed (w c ) and decompressed (to, d ) cloud. Summarising all these corrections then yields the temperature 
during transport 

T = ^ ■ (T im - R c t c ) - R d t d (29) 

Cd c 

with the waiting time t c ( t d ) and heating rate R c ( R d ) in the compressed (decompressed) trap. 
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Chemical potential in the point contact 


Determining the chemical potential using the equation of state implies the assumption of and ideal harmonic trap. 
However the more complex gemoetry at the center of the cloud in the underlying potential of our system due to the 
gate potential and the confinement beams in x— and z —direction must be taken into account. We corrected for two 
effects in estimating the chemical potential in the point contact: 

1. The first correction originates from the fact that the QPC confinement in z expels atoms from the center of 
the cloud, leading to an increased global chemical potential. This confinement is absent during imaging, thus 
it must be corrected for explicitly. The confinement pushes atoms into the reservoirs, increasing the chemical 
potential /x con / with respect to the harmonic trap case p. The correction factor b = p C onf/p is dependant on 
the trapping frequency of the harmonic trap u> and the total atom number N. We estimate b by computing the 
density distribution using the Thomas-Fermi approximation in the full trap based on the homogeneous equation 
of state. For the measurements presented here we find correction factors b between 1.21 and 1.27. 

2. A second correction to the chemical potential arises from finite darkness (~ 99.8%) of the same repulsive beam 
in the QPC, leading to a n additional repulsive optical potential in the nodal line where atoms reside, with an 
amplitude Vfd = 140 nK and a gaussian waist oi Wfd = 30 /im. This is confirmed by a measurement of the 
conductance of non interacting Fermions in the point contact, compared with an ab-initio calculation with the 
Landauer formula. 

The chemical potential and the gap are evaluated in the trap center under disregard of confinement induced effects. 
We are presentaly not aware of calculations of the gap and chemical potential in the presence of 3D scattering and 
weak confinement. 


Numerical derivative of decay curves 


To obtain the current-bias curves shown in the main text, we calculate the numerical derivative of the measured 
AN with respect to time. To this end we extract the local slope using a linear fit over a window of six consecutive 
points. We then shift the window by one point and repeat. Note that the overlap of windows creates correlations 
between adjacent point in the derivative, which leads to oscillating artefacts in the current bias curves. 

From the local slope dAN/dt one obtains the current I: 

_ dN L dN R 1 d(N L - N r ) 1 dAN 

dt dt 2 dt 2 dt [ } 


Exponential fits 


We fit the decay curves at various temperatures (see Figure S3) with an exponential fit of the from 


AN 

IsT 


( t ) 


= n(t) = n^e 


(31) 


with the initial imbalance n o and the decay time r as fit parameters. To test the merit of the exponential fit, we 
deduce the deviation between the fit and experiemental data. For each decay curve only values for t < 3 • r are taken 
into account, and the obtained value \ 2 is divided by the number of points p: 


X 


2 


1 

P 


E 

i,ti < 3 t 


((A N/N)i - n(U)) 2 
—2 


(32) 


Here (A N/N)i is the imbalance form the data at ti, cti is the standard deviation from averaging in AN/Ni and n(t,) 
is the imbalance as obtained from the fit at tj. 
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Current determination 


To obtain the normalized current shown in Figure 4 we compare the initial particle imbalance ANo to the imbalance 
AN after a short time At, determined for each value of the gate potential separately. The mean current / in this 
time window is then given by I = Ajv 'V. Ajv . Normalizing this with respect to the chemical potential bias A/i yields 


I _ Cl 
A ~n ~ AN 


(33) 


which is plotted in Figure 4. Here we made a linear approximation in the compressibility C = AN /A/i, and used the 
harmonic mean AN between AN and ANq to account for averaging over a finite range of the bias. 
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FIG. SI: Tunnelling amplitudes fitted at low temperature as a function of gate potential. For the lowest gate potentials, a 
single tunnelling parameter is fitted for the lowest mode (black rings). For higher gate potentials, three amplitudes are fitted, 
one for the lowest mode and two identical for the next two quasi-degenerate modes (red diamonds). 



FIG. S2: Ratio C/Cq as a function of T/Tp. The veritcal line indicates the critical temperature at fi/kgT = 2.5 



Transport time (s) 


FIG. S3: Decay curves for various temperatures: T = 124 nK (open circle), 135 nK (filled triangle), 145 nK (open triangle), 
170 nK (filled hexagon), 223 nK (empty diamond) and 290 nK (filled square). The solid lines are exponential decay fits (see 
Eq. (311). The data are offset by 0.05 units vertically for clarity. 









